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MIXED-DOMAIN ANALOG/RF SIMULATION 



Field of the Invention 

The present invention generally relates to simulation, and more particularly to 
5 analog (time domain) and RF (time-frequency domain) simulation. 

Background 

The design of an integrated circuit (IC) is an extremely complicated task and 
requires a structured approach. Generally, the IC design flow can be broken down 

10 into the steps of implementation and verification. Implementation usually involves 
the creation of a specification, an RTL model, a transistor-level model, or a gate- 
level netlist. Verification is usually performed by simulators that test the functionality 
of an IC in a software environment before creation of a physical prototype. To 
perform such testing, simulators predict the behavior of a system by using 

15 mathematical models of the system components. Based on the circuit signals, 
different types of simulators are chosen to perform the simulation, such as digital, 
analog, and RF simulators. 

For digital simulation, the modeling begins with a software program that 
describes the behavior or functionality of a circuit. This software program is written 

20 in a hardware description language (HDL) that defines an algorithm to be performed 
with limited implementation details. Designers direct behavioral synthesis tools to 
generate alternate architectures by modifying constraints (such as clock period, 
number and type of data path elements, and desired number of clock cycles). 
Behavioral synthesis tools convert the HDL program into a register transfer level 

25 (RTL) description, which is used for simulation. The RTL description is used to 
ultimately generate a netlist that includes a list of components in the circuit and the 
interconnections between the components. This netlist is used to create the 
physical integrated circuit. Digital simulation usually relates to events where data 
changes on clock cycles and the logical voltage levels are limited to two or three. 
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Analog simulation is used for DC, AC and transient analyses and operates on 
a transistor-level description, which is a full netlist of the physical circuit containing 
detailed architectural and parasitic information. Analog behavioral languages (e.g., 
Verilog-AMS or VHDL-AMS) also exist and allow the user to model an analog 
5 design at a higher level of abstraction. The analog design process includes the 
creation of schematics, layout, and extraction resulting in the netlist used for 
simulation. The analog simulator (e.g., SPICE) uses the netlist as an input to 
formulate non-linear equations for the circuit to solve for the unknown waveform 
V(t). Iterative methods, such as Newton-Raphson, can be used to solve these non- 
10 linear systems of equations. For transient analysis, most analog simulators choose 
time points at which the system is solved in order to meet a desired accuracy 
criteria. Thus the basic difference between analog and digital simulators is that an 
analog simulator considers the voltages and currents in a circuit to be continuous 
quantities, rather than quantized into high/low values, as in digital. This allows 
15 analog simulators to calculate voltage and current as a function of time. 

RF simulators perform both steady-state analysis and modulated steady- 
state analysis (also called envelope analysis). The steady state is defined as the ■ 
state when all initial transients have vanished, and the circuit operates with periodic 
or quasi-periodic large signals. Modulated steady-state analysis is a mixed time- 
20 frequency algorithm dedicated to the simulation of circuits stimulated by (non 

periodic) modulated signals. It efficiently handles the modulation information carried 
by RF signals. A mathematical representation of modulated steady-state signal is 

K 

given by the equation v(t) = ^V K (t)e M{t)t , which includes both time aspects and 

k=-K 

frequency aspects. The output of the modulated steady-state analysis is a time- 
25 varying spectrum. Other analyses are also possible, such as the steady-state small 
signal analysis, steady-state noise analysis, etc. An example use of RF simulators 
is for RF transceivers found in mobile phones. These circuits include low noise 
amplifiers, mixers, filters, oscillators etc., which operate in a RF frequency range 
(e.g., typically from a few 100MHz to 5 GHz). 
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Although a single simulator may have the capability to solve both analog and 
RF, the simulator has two kernels, one for each type of simulation. For example, 
one simulation flow needs to be performed for the analog waveforms and a 
separate simulation flow needs to be performed for the mixed time-frequency 
5 waveforms. Having two simulation flows for separate analog and RF reduces 
efficiency during verification and is difficult when developing high-speed circuits 
where performance is complicated. 

Therefore, it is desirable to provide a mixed-domain simulation system 
capable of handling both analog and RF signal analysis in a single simulation flow. 

10 

Summary 

The present invention provides a simulation environment wherein both 
analog and RF signals are simulated in a single flow by a mixed-domain simulator. 

15 The simulator includes a simulator kernel with an analog solver and an RF solver to 
allow both analog- and RF-type of signals to be solved in an interrelated fashion. 
The simulator may also include a partitioner that divides the circuit into various RF 
and analog modules to be solved. User input may control the partitioning process, 
but the simulator may refine the partitions and/or generate sub-partitions to improve 

20 performance and provide a higher probability of convergence. 

These features and others of the described embodiments will be more readily 
apparent from the following detailed description, which proceeds with reference to 
the accompanying drawings. 

25 Brief Description of the drawings 

FIG. 1 shows a mixed-domain simulator system used to simulate circuits 
including both analog and RF signals. 

FIG. 2 is a flowchart of a method for simulating mixed-domain circuits 
according to the invention. 
30 FIG. 3 is a detailed flowchart of a method for simulating mixed-domain 

circuits. 
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FIG. 4 is a detailed flowchart of a method for solving RF and analog partitions 
of FIG. 3. 

FIG. 5 shows an example circuit partitioned into analog and RF sections. 
FIG. 6 provides a flowchart of a method for partitioning. 
5 FIG. 7 shows a network environment in which the system may be used. 

FIG. 8 shows a flowchart of a method for simulating over the network of FIG. 

7. 

Detailed Description 

10 

FIG. 1 shows a mixed-domain analog/RF simulator environment 10 including 
a simulator 12. The simulator 12 has a single simulator kernel 14 able to simulate 
both analog and RF signals in a single simulation flow. This concept is shown 
diagrammatically by including within the simulator kernel 14 separate analog and 
15 RF solvers 16, 18, respectively. The analog solver 16 formulates non-linear 

equations for the circuit to solve for an unknown waveform V(t), where V(t) is a time- 
domain signal. RF solver 18 is used to perform modulated steady-state analysis of 
RF signals. For modulated steady-state analysis, the RF solver might represent the 

K 

RF signal using the following formula: v(0 = ^V K {t)e j(0k{t)t , that includes both time 

k=-K 

20 aspects and frequency aspects. The simulator may also include a partitioner 19 to 
divide the circuit into modules that are separately solved to make the system of non- 
linear equations manageable and, thereby, allowing faster simulation time. A 
relaxation tester 21 performs one-step relaxation to test convergence of the solved 
partitions. Any desired relaxation method may be used, such as that disclosed in 

25 US Patent No. 4,791,593 to Bernard Hennion, which is hereby incorporated by 
reference. 

A circuit description that includes components used with both analog and RF 
signals is provided to the simulator 12 through a netlist 20, which may take many 
forms. For example, the netlist may be a transistor-level description, which is a full 
30 netlist of the physical circuit containing detailed architectural and parasitic 
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information or the netlist may be derived from an analog behavioral language (e.g., 
Verilog-AMS or VHDL-AMS). In any event, the netlist normally contains a list of 
components in the circuit, the parameter values, and connection node names. The 
simulator 12 also receives control statements 22 directly in the netlist or through a 
5 user interface (not shown). The control statements may be used in a variety of 
ways, such as allowing the user to control the type of simulation to be performed, 
the stimuli to be used, which signals are displayed, partitioning, etc. 

The simulator 12 uses the netlist 20, the control statements 22, and one or 
more databases, such as an RF library 24 and analog library 26, to simulate the 

10 circuit and provide simulation results 28 in any desired form, such as on a display, 
in an output file, etc. Although two databases 24, 26 are shown, they may be 
combined or further subdivided. 

FIG. 2 is a flowchart of a method for performing the simulation using the 
simulator 12 of FIG. 1. In process block 40, the simulator 12 reads the netlist that 

15 includes a circuit description containing both components used with analog signals 
and components used with RF signals. In process block 42, the simulator 12 
simulates the circuit using a single-flow simulation that is used to verify the correct 
operation of the circuit. The single-flow simulation is performed by a single 
simulator kernel 14 (FIG. 1) and as further described below. In process block 44, 

20 the simulation results are outputted in any desired manner. For example, the 
simulation results of the RF signals may be displayed at the same time as the 
simulation results of the analog signals, wherein both the RF and analog signals 
were obtained during the same simulation flow. The single simulation flow allows 
for an interrelated solution where solutions of analog partitions affect the analysis of 

25 RF partitions and vice versa. 

FIG. 3 shows a more detailed flowchart of process block 42 (Fig. 2) for 
simulating a mixed-domain analog and RF circuit using the simulator 12. In process 
block 46, the circuit is partitioned into analog and RF modules. Further details of 
the partitioning process are described below in relation to FIGs 5 and 6. In process 

30 block 48, the time variable T is initialized to zero. The variable T stores the current 
time being evaluated because the simulator 12 simulates based on variable 
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increments of time H. Thus, process block 48 initializes T to start at the beginning. 
The single simulation flow then continues with a loop that starts at decision block 
50, where the variable T is analyzed to determine if a stopping point of TFINAL is 
reached (TFINAL being the maximum time to be simulated as defined by the user). 
5 If so, the simulation is complete as shown at process block 52. Otherwise, the 
simulation continues at process block 54 where the time increment H is defined and 
a solution predicted. For example, a solution is extrapolated based on the chosen 
time increment H and one or more of the previous time simulation passes. If a 
resolution failed to converge, a smaller time increment can be used to increase the 

10 probability of convergence. At process block 56, the current time T is incremented 
by H. At process block 58, the RF and analog partitions are solved as described 
further below in relation to FIG. 4. At process block 60, a check is made for 
convergence. If there is convergence, the loop is continued as shown by arrow 62. 
Otherwise, the current time T is decremented at process block 64 and the process 

15 continues as shown by arrow 66. In the case of non-convergence, process block 54 
will take into account that convergence was not successful and will modify H 
accordingly to increase the probability of convergence. 

FIG. 4 shows further details of process block 58. At process block 70, a loop 
index i is set to one. In process block 72, an RF or an analog partition designated 

20 by index i has its boundary nodes positioned. The boundary node positioning 

specifies the boundary node values as fixed values or spectrums to the partition so 
that the partition can be solved using Newton-Raphson techniques. In the case of a 
node that is an analog node coupled to either an analog partition or an RF partition, 
the boundary node is set to a fixed value. In the case of a node that is an RF node 

25 coupled to an RF partition, the boundary node is set to a fixed spectrum defined by 
the V K (k=-K to K). In the case of a node that is an RF node coupled to an analog 
partition, the boundary node may be set in at least two ways: In a first fast mode it 

K 

is set to V 0 (t) (obtained by limiting K to 0 in the summation v(t) = ^V K {t)e ja>k{t)t ). In 

k=-K 

K 

a more accurate mode, it is set to v(t) ( v(t) = ]T V K (r)e M(0 ' , computed over the 
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entire summation). In process block 74, with the boundary nodes held fixed, the 
partitions are solved using Newton-Raphson or other well-known techniques. If the 
partition is an RF partition, the simulator 12 uses the RF solver portion 18 of the 
simulator kernel 14 to solve the system of equations defined by F(V K ) = 0, as is well 
5 known in the art. Alternatively, if the partition is an analog partition, the analog 
solver portion 16 of the simulator kernel 14 is used to solve the system of equations 
defined by F(v)=0. In any event, there is an interrelationship between solving of the 
RF partitions and the analog partitions through the boundary nodes and boundary 
components. 

10 At decision block 76, divergence is checked. If there is divergence, then the 

loop exits at 78 and the process continues at decision block 60 of FIG. 3. Thus, 
returning briefly to FIG. 3, convergence at 60 should also be answered in the 
negative when there is divergence at 76 in FIG. 4. Therefore, the time is 
decremented (block 64), and a new (shorter) time increment is chosen at block 54 

15 to provide a higher probability of convergence. 

Returning to FIG. 4, if there was not divergence at decision block 76, then a 
check is made at decision block 80 whether all of the partitions have been solved. If 
no, then index i is incremented (process block 82) and the loop continues as shown 
by arrow 84. If in decision block 80, all of the partitions have been solved, then 

20 overall convergence is checked by relaxation tester 21 . Decision block 86 

determines whether there was convergence and, if yes, then the routine is exited at 
88. If there is no convergence, then the process is repeated as shown by arrow 90. 

FIG. 5 is an example circuit 100 divided into partitions shown generally at 
102. Although only a small number of partitions are shown, there may be any 

25 number of partitions based on the circuit. The partitions 102 are divided into two 
general sets of nodes and components: RF 104 and analog 106. Boundary nodes 
(shown by dots such as shown at 108) and boundary elements 1 10 are shown 
between the partitions 102. Thus, a partition is a set of nodes and a set of 
components, wherein any given node belongs to only one partition. A boundary 

30 node is a node connected to a (boundary) component that is connected to at least 
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one node in another partition. A boundary component is a device connected to 
boundary nodes and is shared by at least two partitions. 

FIG. 6 shows a flowchart expanding on process block 46 of FIG. 3. In 
process block 120, user input is received regarding how to partition the circuit. The 
5 user can provide control statements 22 about how the circuit should be partitioned 
generally. However, the simulator refines this information to create the actual 
partitions, such as which RF nodes are placed in which RF partition, etc. In 
process block 122, the simulator partitions the analog and RF components and 
nodes based on the user input. In process block 126, the partitioner 19 examines 

10 and automatically refines the partitions (and nodes) to encourage convergence. 
The simulator 12 uses information in the netlist, the user-provided information, and 
its own refinements to form a partitioned circuit 100, such as shown in FIG. 5. In 
process block 128, the simulator 12 may create further sub-partitions to speed-up 
simulation. Finally, in process block 130, the partitions 102 are ordered to define 

15 the resolution flow. Generally, the partitions are ordered to follow signal flow 

through the circuit. For example, a partition coupled to an input signal may be first 
in order. In any event, the partition order corresponds to the index number i of FIG. 
4. 

FIG. 7 shows that portions of the system 10 may be applied to a distributed 
20 network, such as the Internet. For example, a server computer 150 may have an 
associated database 152 (internal or external to the server computer). The server 
computer is coupled to a network shown generally at 154. One or more client 
computers, such as those shown at 156 and 158, are coupled to the network to 
interface with the server computer using a network protocol. 
25 FIG. 8 shows a flow diagram using the method on the network of FIG. 7. In 

process block 170, the netlist is sent from a client computer, such as 156, to the 
server computer 150. In process block 172, the circuit described in the netlist is 
simulated using a single-flow, mixed-domain simulator 12. In process block 174, 
the results are returned to the client computer via the network 154. In process block 
30 1 76, the client computer displays or otherwise provides the simulation results to the 
user. 
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Having illustrated and described the principles of the illustrated 
embodiments, it will be apparent to those skilled in the art that the embodiments can 
be modified in arrangement and detail without departing from such principles. 

Although boundary components are shown as being connected to only two 
5 partitions, they may be connected to any desired number of partitions. 

Although it is described that the user defines some information to control 
partitioning, the partitioning may be fully automated by the simulator. For example, 
the simulator can analyze whether an input stimulus is RF or analog and control 
partitioning based on that analysis. Additionally, the signal types at the boundary 
10 nodes can be used to define the partitioning. 

In view of the many possible embodiments, it will be recognized that the 
illustrated embodiments include only examples of the invention and should not be 
taken as a limitation on the scope of the invention. Rather, the invention is defined 
by the following claims. We therefore claim as the invention all such embodiments 
15 that come within the scope of these claims. 



